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The inclusive B — > X s £ + £~ decay rate in the large q 2 region (q 2 > miL) receives significant 
nonperturbative corrections. The resulting uncertainties can be drastically reduced by normalizing 
the rate to the B — > X u £u rate with the same q 2 cut, which allows for much improved tests of 
short distance physics. We calculate this ratio, including the order nonperturbative correc- 

tions and the analytically known NNLO perturbative corrections. Since in the large q 2 region an 
inclusive measurement may be feasible via a sum over exclusive states, our results could be useful 
for measurements at LHCb and possibly for studies of B — > X d £ + £~ . 



I. INTRODUCTION 

The b —> s £ + £~ process plays an important role in 
making overconstraining measurements of CKM matrix 
elements and searching for physics beyond the Stan- 
dard Model (SM). This decay has been observed both 
in inclusive B — > X s £ + £~ [jl], ^ and exclusive B — > 
K^h+fr |, H transitions. The inclusive B -> X s i+l~ 
decay rate can be calculated in a systematic expansion if 
one ignores the J/ip and ip' resonances. It has thus been 
advocated to compare calculations and measurements of 
the (differential) rate for q 2 < m 2 /^, and q 2 > m^, , which 
we shall refer to as the small q 2 and large q 2 regions, 
respectively. Here q 2 — (p t + + Pi-) 2 is the dilepton in- 
variant mass, and in practice the q 2 regions are chosen 
as q 2 < 6 GeV 2 and q 2 > 14 GeV 2 . 

The measurements in the two regions are complemen- 
tary, as they have different sensitivities to short distance 
physics, the main theoretical uncertainties have different 
origins, and the experimental challenges are also distinct. 
The most important operators for B — > X s £ + £~ are 
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that the B — > X s £ + £~ rate in the large q 2 region has a 
large hadronic uncertainty [Q, ||, The reason is that 
the operator product expansion becomes an expansion in 
Aqcd/(™6 - \fq*) instead of A QC d/™&- 

Our main point is that this uncertainty can be dras- 
tically reduced by comparing measurements and calcula- 
tions of the ratio 
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with the same lower cut q 2 > q^ in the b — > s and b — > u 
decays. 1 The nonperturbative corrections related to the 
dominant Og and Oiq contributions are the same as for 
the semileptonic rate. Thus, as explained below, nonper- 
turbative effects in the ratio in Eq. (Q) are suppressed 
near maximal q 2 by 
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which is nearly an order of magnitude. The scheme we 
use for the Wilson coefficients C7 9 10 (L2j will be defined 



in Sec. III. Their SM values are 



The operator Or is important in B — > X s £ + £~ at small 
q 2 due to the 1 jq 2 pole from the photon propagator (and 
it dominates the B — > X s j rate). At large q 2 , however, 
the O7 contribution is small. Compared to small q 2 , the 
rate in the large q 2 region has a smaller renormalization 
scale dependence and m c dependence ||. Although the 
rate is smaller at large q 2 , the experimental efficiency is 
better (l], |^|. Moreover, requiring large q 2 constrains the 
X s to have small invariant mass, mx 3l which suppresses 
the background from B — ■> X c £~v — > X s £ + £~vv. To 
suppress this background at small q 2 , an upper cut on 
mx 3 is required, complicating the theoretical description 
due to the dependence of the measured rate on the shape 
function ||, which is absent at large q 2 @, §. 

Despite these advantages, the large q 2 region has been 
considered less favored. The 1/mf corrections are not 
much smaller than the 1 / m 2 , ones M , so it is often stated 



C 9 = 4.207, Ci = -4.175, C 7 = -0.2611 



(4) 



We calculate in this paper the ratio TZ(qg), which al- 
lows one to translate the measured B — > X s l + l~ and 
B — * X u £v rates in the large q 2 region to a precision 
constraint on the Wilson coefficients times the CKM el- 
ements in Eq. (|2|). The normalization of 72.((?q) is chosen 
such that at lowest order and in the limit IC7 /Cg^io | <C 1, 
Tl{ql) = Cf + Cf a . Hereafter we assume the SM and 
neglect the strange quark and lepton masses. 



1 This was noted as a remote possibility in Ref. |8l, but was sub- 
sequently forgotten even by those authors. Tne experimental 
prospects have improved sufficiently that such a study may be 
possible in the near future. 
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II. THE q 2 SPECTRA TO ORDER 

The nonperturbative corrections to the q 2 spec- 
trum are calculable in an operator product expan- 
sion (OPE) [Oh The first corrections appear at 
0(Aq CD /m 2 ) |{I4t P^| - They are parameterized by 
two nonperturbative matrix elements, Ai and A2. At 



0(Aq CD /?tt|) there are two new local matrix elements, 
Pi and p 2 , four time-ordered products, 7i_ 4 jl6|, |l7| , and 
process dependent matrix elements of four-quark opera- 
tors, fi @|1]||2(J. 

The q 2 spectrum up to 0(A^ CD /m 3 ) for B — » X u lv is 
given by || H 
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where s = q 2 /m 2 , and 1/(1 - x)+ = lim £ ^ o [0(l - z - e)/(l - x) + S(l - x - e) lne]. For B -> X s t+l~ [g, [TJ || 0, 
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The nonperturbative parameters in Eqs. (0) and (||) are 
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For our purposes, the % can be absorbed into Ai >2 . In the 
total rate and the q 2 spectrum, Ai enters proportional to 
the b quark decay rate, and the p 2 contribution is pro- 
portional to A 2 . Hence, the important nonperturbative 
parameters for the q 2 spectrum are A 2 , p\, and f U}S . 

The value of A 2 is known fairly precisely, A 2 = (m 2 Bt — 
m 2 3 )/4 ~ 0.12 GeV 2 . To estimate pi, the equations 
of motion can be used to relate the relevant operator 
to a four-quark operator. Using the vacuum saturation 
model, pi = (27ra s /9)/|m B ~ (0.4 GeV) 3 @ H- The 
fits to the B — > X r lv shape variables are sensitive to p\ 
and prefer a larger central value J2^, |24| with significant 
uncertainties. We shall use p\ — (0.1 ± 0.1) GeV 3 . 

The four-quark operator contributions, f u and f s 
(sometimes called weak annihilation, though the light 
quark flavor need not match the flavor of the specta- 
tor quark), depend on the final state and on the flavor 



of the decaying B meson. They contribute near maximal 
q 2 , and their contribution has only been derived for the 
total rate fl8j pl| p0[ |2l| and the lepton energy spec- 
trum [fl8[ |26| . However, the four-quark operators have 
to be consistently included in the OPE for the fully dif- 
ferential spectrum. This affects the matching for pi in a 
nontrivial way at s = 1, replacing the singular pi/(l — s) 
terms present in the earlier literature by the plus distri- 
butions in the last lines of Eqs. ([5]) and (|^). Apart from 
this unambiguous regularization of the singular integrals 
at s = 1, our result in Eq. (|^) agrees with Ref. ||. 

The values of f u and f s are poorly known. They are 
important, since they are enhanced by a loop factor, 
16tt 2 . In the notation of Ref. (2fJ, f u = 2n 2 j 2 B m B (B 1 - 
B 2 ), where B\ 2 are phenomcnological "bag parameters". 
In the vacuum saturation model, B\ — B 2 — 1 in charged 
B decay and B\ — B 2 = in neutral B decay. This gives 
a significant suppression with large uncertainty, since 
the accuracy of the model is poorly known. Because of 
this sensitivity to cancellations between nonperturbative 
quantities with comparable magnitudes, the estimates of 
f u are uncertain. On general grounds one expects in 
charged B decay to be greater in magnitude than /° 
in neutral B decay, since in the former case the spectator 
flavor matches the flavor of the light quark in the four- 
quark operator. The assumption \Bi — B 2 \ = 0.1 (2TJ] for 
charged B decay leads to |/^| = 0.4 GeV 3 . This would 
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FIG. 1: Left: The impact of the corrections proportional to A2 (suppressions) and pi (enhancements) on the Cg + Ci , C7C9, 
and Cf contributions to the partonic B — * X s £ + £~ rate integrated over so < s < 1. Right: The partonic B — ► X s £ + £~ rate for 
so < s < 1 divided by the Cg + C10 contribution (the latter is proportional to the B — » X u £v rate), at lowest order (dotted), 
including the A2 terms (dashed), and including both A2 and pi corrections (solid). (Note the different scales.) 



give a large uncertainty in 7Z(qq), and we discuss next 
how it can be reduced. 

Flavor ST/ (3) symmetry implies /° « /°. One also ex- 
pects fg w , though this requires assumptions beyond 
ST/ (3) (as there are two singlets in 3 x 3 x 3 x 3). Thus, 
one may use in Eq. @ the average of the charged and 
neutral B — > X s £ + ^~ rate and the neutral £?° — » X u tv 
rate in the presence of the same g 2 cuts. 2 Separately 
measuring B° — > AT„£p and i? 1 * 1 — > in the large 

q 2 region is important, even without the additional rea- 
sons discussed here, for the determination of \V u b\ and 
to constrain weak annihilation. We anticipate that this 
separation will be available by the time the B — ► X s £ + £~ 
rate in the large q 2 region is precisely measured. 

To illustrate the slow convergence of the OPE at large 
q 2 , we show in the left plot in Fig. [j] the effect of the 
dominant 1/w 2 and corrections, proportional to 

A2 and pi, respectively. For each of the Cf + C 2 , CrCg, 
and Cj contributions, we plot the free quark decay rate 
plus either the A2 (curves below unity) or the p\ (curves 
above unity) terms integrated over so < s < 1, normal- 
ized to the free quark decay rate. The A2 term, which 
is about a —2% correction to the total rate, makes the 
spectrum negative for s > 0.9, and the integrated rate 
negative for s > 0.8. The terms proportional to p\ give 
a correction of comparable magnitude and opposite sign 
(taking p\ = 0.1 GeV 3 ). Even at the rather low cut, 
sq = 0.6, the correction to the rate proportional to A2 
is about —21%, while that proportional to pi is about 
+17%. Moreover, these estimates do not include the /j 
terms discussed above. Thus, even for the rate integrated 



2 For B - X d t+l-, 5(7(3) implies f° « /± and /± « f°. If 
semileptonic B a — * Xv,lv decays can be measured in the large q 2 
region, that would open up yet another set of applications. 



over the entire large q 2 region, it is not clear how to assign 
a robust uncertainty. 

The important point is that while these corrections 
are large, the terms proportional to C| + Cf dominate 
at each order, and they are identical to those that occur 
in B — » X u £u. We present a detailed numerical study in 
the next section. To illustrate the point, the right plot in 
Fig. [l] shows the partonic B — > X s £ + £~ rate integrated 
over so < s < 1 divided by the C| + C 2 contribution, at 
lowest order (dotted), including the A2 terms (dashed), 
and including both A2 and p\ corrections (solid). The 
blowup of the dashed curve near sq ~ 0.8 is unphysical 
and occurs because the rate truncated at C(Aq CD /to 2 ) 
becomes negative at slightly different values of sq for the 
C|+C 2 , C7 Cg, and C 2 contributions, as can be seen in the 
left plot. As expected, for so — ► 1 the ratio is 0.88 at each 
order, i.e., unity minus the 0.12 correction in Eq. (0). 



III. RESULTS FOR K(q%) 

To organize the different short-distance contributions, 
and combine the nonperturbative and the next-to-next- 
to-leading order (NNLO) perturbative corrections, we 
use the scheme introduced in Ref. p3|. The Wilson co- 
efficients £7,9.10 are defined as 

C 7 = C 7 {p) [m b (p)/m b ] + ... , 
C 9 = C 9 (p) + ... , 

Cio = C w . (9) 

where the ellipses denote a minimal set of perturbative 
corrections, such that are fi independent and real 
in the SM (which is automatic for Cm). We use the IS 
scheme |27j , which improves the behavior of the pertur- 
bative expansions (for the semileptonic q 2 spectrum both 
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the a 2 (3g || and full a 2 (23] corrections are known). Con- 
sequently, we use nib = everywhere, except for the 
MS 6-quark mass, m&(/z), in Eq. ([)[), which is renormal- 
ized together with CV(/i). 

The inclusive decay rate is expressed in terms of the 
effective Wilson coefficients lOl 
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which are defined such that all terms on the right-hand 
side are separately u, independent to the order we are 
working at. We view the coefficients C 7 $.\q as the param- 
eters sensitive to physics beyond the SM, which should 
be extracted from experimental data and compared with 
their SM predictions. Even if they receive significant new 
physics contributions, the functions F 7 ^(q 2 ) and G 7 ,g{q 2 ) 
are likely to be dominated by the SM. 

The ^7,9 (g 2 ) terms in Eq. ( |l0| ) contain contributions 
from the remaining Oi_6.8 operators in the effective 
Hamiltonian, for which we employ a partial NNLO treat- 
ment. We use the Wilson coefficients at 0(a s ) [^9| 
pl| , |32| , j33j i _but only keep the O(a s ) contributions to 
Fg^q 1 ) [fljUl. (Note that F 7 (q 2 ) vanishes at order a° s .) 
We cannot include the 0(a s ) corrections to Fr ! g(q 2 ), be- 
cause the dominant O1.2 contributions are only known 
analytically in the small q 2 region |35| |. They have been 
evaluated numerically for large q 2 , and lead to a reduced 
scale dependence and central value || . 

The G 7t g(q 2 ) terms in Eq. (|l(J) contain the calcu- 
lable Aq CD /m^ nonperturbative corrections associated 
with intermediate cc loops f36| , |37f . They can be in- 
cluded this way for any differential rate |l2| ], provided 
0(a s /m 2 , f/m^) cross terms are neglected. They affect 
Tl(q 2 ) slightly below the 1% level in the large q 2 region. 

Thus, the q 2 spectra, including up to NNLO perturba- 
tive and l/m^ nonperturbative corrections, are 
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The power suppressed corrections, dT}J™ / dq 2 , are given 
by the terms in Eqs. (g) and (§) proportional to A12, 
pi t 2, and f u>s , appropriately replacing C 7 ,g by C\ n f as in 
Eq. (P). The functions Sl ij (s), 



K 9 ( S ) + 2 S 4 9 ( S )], 
[ S u 7 L 7 ( S ) + 2to 7 T 7 ( S )], 
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parameter 


central value 


uncertainty 


H [GeV] 


4.7 


+4.7 
-2.35 


m b [GeV] 


4.7 


±0.04 


m c [GeV] 


1.41 


±0.05 


A 2 [GeV 2 ] 


0.12 


±0.02 


„ 1"/^ . 7"3l 

pi [GeV J 


0.1 


±0.1 


At [GeV 3 ] 





±0.4 


f° - fs [GeV 3 ] 





±0.04 


/2 + fs [GeV 3 ] 





±0.2 



TABLE I: Central values and ranges of input parameters. The 
parameters Ai, pi, % are irrelevant for this work and are set 
to their central values, Ai = —0.27 GeV 2 and pi = % — 0. 



contain the 0(a s ) corrections to the matrix elements of 
the OiOj contribution fl], |J, H H Hi converted to the 
IS scheme, with w^ T (s) given in Ref. jl2). We neglect 
finite bremsstrahlung corrections, r brcms , associated with 
01-6,8, because they are negligible at large q 2 |4jJ . 

The perturbative uncertainty due to the choice of 
renormalization scale for a cm , which appears in the pref- 
actor of K(ql) in Eq. (§), can be eliminated by includ- 
ing the relevant higher order electroweak corrections to 
the B — > X s £ + £~ rate, which have been studied in 
Refs. [||, ||] in the small q 2 region. 

We present our result for the SM prediction for two 
different values of g 2 ,, as it is an open question what its 
most suitable choice is. Using the input values in Table Q, 
with all other inputs as in Ref. p2[, we obtain 



K(U GeV 2 ) 
^(15 GeV 2 ) 



+ 4 
± 1.06 C 9 
£9 4 C 2 4 
±0.97 C 9 



4.79 C 2 + 4.3IC7C9 
-2.24C 7 4 0.95, 
4.27Cf + 4.10C 7 C 9 
- 1.91 C 7 4 0.93. (13) 



Using the central values of the C; in Eq. ([!]) and evalu- 
ating the uncertainties by varying the parameters within 
their ranges given in Table [jj we find 

^(14 GeV 2 ) = 35.55 (l ± 0.046 [/u s] ± 0.012 [A2iPl] 
± 0.054 M ±0.030 [Ci] ), 

K(15 GeV 2 ) = 35.42 (l ± 0.065 [/u s] ± 0.016 [Wl] 

± 0.051 M ±0.030 [Ci] ). (14) 

Eqs. (|lf) and (@) can be directly compared with experi- 
mental measurements to constrain the Wilson coefficients 
(mainly C| + C 2 ) and test the standard model. 

The first two uncertainties in Eqs. ( |l4| ) are due to 
nonperturbative corrections. They are shown in Fig. |2, 
where we plot 7£(gg)/7£(14 GeV 2 ) as a function of q$. 
The most important uncertainty is due to the four-quark 
operators (weak annihilation), in particular, the differ- 
ence f u — f s , which does not cancel in the ratio 7£(<7q). 
The green (wide light) region in Fig. corresponds to 
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FIG. 2: Nonperturbative uncertainties in H{ql) /TL{1A GeV 2 ). 
The black curve is the central value. The green (wide light), 
blue (dark), and orange (narrow light) shaded bands show the 
uncertainties from f u ,s, pi, and A2, respectively. 



FIG. 3: Perturbative uncertainties in 7£(gnY/7£(14 GeV 2 ). 
The black curve is the central value (as in Fig. H) . The orange 
(narrow light) and green (wide light) bands show the p de- 
pendence (see text). The blue (dark) band is the uncertainty 
from £7,9,10. 



allowing l/o - f„\ < 0.04 GeV 3 and |/° + f s \ < 0.2 GeV 3 , 
which is appropriate if the average of the B° and B^ rare 
decay data is compared with the B° semileptonic data. 
Fixing = which holds in the SU(3) limit, and 
varying them together in the range |/° s | < 0.1 GeV 3 , 
gives a very small uncertainty, below 1%. If B° — > X u £v 
is not measured separately from B ± — > X u £v, the un- 
certainty in 7\L(14GeV 2 ) increases to above 20% (taking 
|/ u | < 0.2 GeV 3 and f s — as an estimate for the average 
of the charged and neutral modes). 

Even with our conservative range for p\ , the combined 
uncertainty from p\ and A2 is very small, at the 1% level. 
The allowed variation of A2 accounts for T3 4, which it ab- 
sorbs in Eq. (|J), and for p 2 , which it is proportional to 
in Eqs. (JsJ) and (^|). The individual uncertainties from 
pi and A2 are shown in Fig. ^ as the blue (dark) and or- 
ange (narrow light) bands, respectively. The total uncer- 
tainty from nonperturbative corrections, estimated con- 
servatively, is only about 6%. Without normalizing to 
the B —> X u lv rate with the same cut, the uncertainty 
in the B -» X s £+£- rate for q 2 > 14 GeV 2 due to A 2 , f s , 
and pi would be 4%, 9%, and 21%, respectively. 



The remaining two uncertainties in Eqs. (14) are due to 
the perturbative corrections, and are illustrated in Fig. |^. 
The variation in p,, shown by the green (wide light) band, 
is almost entirely due to the 0(oP 8 ) pieces in F 9 (q 2 ), be- 
cause their p dependence does not cancel in lZ(q 2 ). We 
expect that it will be reduced to below 1% by includ- 
ing the full NNLO result. This is shown by the narrow 
orange (narrow light) region in Fig. |^, obtained by includ- 
ing the subset of a s corrections to Fg(q 2 ) which cancels 
its leading p dependence (see Eq. (A5) in Ref. |l2"|). Fi- 
nally, the blue (dark) band shows the uncertainty from 
£7,9,10, which includes their residual p, dependence and 
dependence on the electroweak matching scale (mainly 
affecting Cg), as well as their dependence on the top- 
quark mass (mainly affecting Cio). The uncertainty in 
Fj : g(q 2 ) due to the top quark mass is negligible and that 



due to the electroweak matching scale is well below 1%. 

The uncertainties due to other input parameters are 
much smaller than those shown in Eq. (|14|). The uncer- 
tainty from m c is well below 1%, about the size of that 
of A2, shown by the orange (narrow light) band in Fig. 0. 
The TOfc dependence is negligible, because it almost com- 
pletely cancels between the numerator and denominator 
of 7\L(g 2 ). This is another significant advantage of consid- 
ering the ratio 7?.(g 2 ). Both integrated rates with a cut 
at q 2 — 14 GeV 2 scale roughly as ml 3 , yielding a 11% 
uncertainty. Normalizing the B — > X s l + l~ rate to the 
total B — > X u £v rate (proportional to mf), would still 
leave about an m\ dependence (and even stronger if nor- 
malized to B — * X c £v), which would give an additional 
7% uncertainty compared to our results. 



IV. CONCLUDING REMARKS 

In this paper we pointed out that the theoretical un- 
certainty of the B — » X s £ + £~ rate in the large q 2 re- 
gion, which is dominated by nonperturbative uncertain- 
ties, can be significantly reduced by normalizing the 
B — ► X s l + £~ rate to the B — > X u £v rate with the same 
q 2 cut. Fully exploiting this proposal requires the exper- 
imental separation of B° — > X u £v and B ± — * X u £v in 
the large q 2 region to eliminate the largest part of the 
four-quark operator (weak annihilation) contributions. 

In much of the theory literature B — ► X s £ + £~ has been 
normalized to B — > X c lv and sometimes to B — > X u £v. 
However, to achieve the most reduction of theoretical 
uncertainties, related to both the matrix elements of 
the higher dimension operators in the OPE and the 
value of rrib, the normalization should be done using the 
B — > X u £v rate with the same cuts as in B — > X s £ + £~ . 
This holds both in the small q 2 region || , and especially 
in the large q 2 region studied in this paper. 



6 



An uncertainty we have not addressed explicitly is due 
to higher cc resonances, ^(35) and above, but their con- 
tributions are much smaller than those of the if) and ip' 
and are not expected to introduce a significant uncer- 
tainty. Improvements in both theory and experiment 
will determine what is the optimal choice of the lower 
limit, °f the large q 2 region to minimize the to- 
tal uncertainty. If g§ is increased much above rr&,, 
one should get concerned about quark-hadron duality. 
Even for q 2 = 14GeV 2 the range of hadronic invari- 
ant masses summed over is only mx„ < 1.53 GeV (and 
m Xs < 1.41 GeV for q 2 = 15 GeV 2 ). Hadronic r decay 
data [Q indicate that duality may already be a good 
approximation at these values. 

In the small q 2 region, in addition to the q 2 spectrum 
and the forward-backward asymmetry, dApB/dq 2 , one 
can constrain a third linear combination of the Wilson 
coefficients by splitting the rate into "transverse" and 
"longitudinal" components At large q 2 this is not 

advantageous, since both are dominated by the Cf + C 2 
contribution, and thus yield very similar constraints. 
Our results can also be applied to the forward-backward 
asymmetry, Afb, which at large q 2 is mainly sensitive 
to CgCio- As noted in Ref. M, the OPE at large q 2 ap- 
pears to be significantly better behaved for Afb than 
for the q 2 spectrum. However, if only the normalized 
forward-backward asymmetry, (dApB/dg 2 )/(dr/d<2 2 ), is 
measured in the large q 2 region (which may be the least 
difficult experimentally), it would inherit the large un- 
certainty of the rate. In this case, one could normalize 
this measurement to the rate or the normalized Afb in 
B — > X u tv with the same q 2 cut. The latter should be 
accessible with the large samples of fully reconstructed 



B decays used to extract \V u b\- 

It is not yet known if an inclusive study of B — > 
Xd£ + l~ can ever be carried out, but it may be less dif- 
ficult in the large q 2 than in the small q 2 region. At 
large q 2 a semi-inclusive experimental analysis might be- 
come feasible at a super B factory or even at LHCb. The 
methods discussed in this paper are clearly applicable to 
this decay as well. Moreover, in the large q 2 region the 
exclusive rates can be understood model independently 
using continuum methods |4f| or lattice QCD. 

To summarize, we showed that it is possible to gain 
theoretically clean short distance information from the 
large q 2 region of B — > X s £ + i~. For this, the experi- 
mentally most important input, in addition to precisely 
measuring B —> X s l + £~ at large q 2 , is the measurement 
of B a — > X u lv with the same q 2 cut without averaging 
with B ± — > X u lv. On the theory side, it is desirable 
to include the full NNLO calculation, F 7 ^{q 2 ) at 0(a s ), 
which will largely reduce the perturbative uncertainties 
in Fig. [| and leave us with a prediction for 7£(g 2 ) in the 
14-15 GeV region with a theoretical uncertainty about 
the 5% level. 
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